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The time-convolutionless master equation for the electronic populations is derived for a generic 
electron-phonon Hamiltonian. The equation can be used in the regimes where the golden rule 
approach is not applicable. The equation is applied to study the electronic relaxation in several 
models with the finite number normal modes. For such mesoscopic systems the relaxation behavior 
differs substantially from the simple exponential relaxation. In particular, the equation shows the 
appearance of the recurrence phenomena on a time-scale determined by the slowest mode of the 
system. The formal results are quite general and can be used for a wide range of physical systems. 
Numerical results are presented for a two level system coupled to a Ohmic and super-Ohmic baths, 
as well as for a a model of charge-transfer dynamics between semiconducting organic polymers. 



I. INTRODUCTION 



II. HAMILTONIAN AND THE INITIAL STATE 



The Pauli master equation and the Redfield equation 
have long been applied to the study of relaxation dynam- 
ics in electronic systems 0, 0, Q . These equation were 
originally derived using semi- heuristic arguments. Hww- 
ever, with the advent of projection operator techniques 
m 1960's iii it became clear that such master equa- 
tions could be obtained from formally exact generalized 
master equations through a series of well defined approxi- 
mations. Such approximations are usually justified when 
the system of interest is coupled to a thermostat with a 
large number of degrees of freedom. In this limit, the re- 
sultant transition probabilities in the Pauli master equa- 
tion are identical to those obtained by the less rigorous 
golden rule approach . 



A wide class of the electron-phonon systems can be 
described by the following Hamiltonian {h = 1) 



H 



e„|n) (n| -I- 



\n) (m|(aj -|- a^) + } ^ujtalai.jl) 



Here |n)'s denote electronic states with vertical energies 
e„, a| and Ui are the creation and annihilation operators 
for the normal mode i with frequency tOi , and gnmi are the 
coupling parameters of the electron-phonon interaction 
which we take to be linear in the phonon normal mode 
displacement coordinate. 

We can separate H into a part that is diagonal with 
respect to the electronic degrees of freedom, 



|7i)(7i| gnni\n){n\{al + Ui) -|- Wjajai,(2) 



and an off-diagonal part V 

V = 5mm I") ("^1 (4 + Oi): 



(3) 



The case when the thermostat is mesoscopic (i.e. with 
a finite number of modes) is less well studied. In general, 
for this case the time independent transition probabilities 
of the golden rule type may not be well defined since the 
kernels of the generalized master equations do not vanish 
in the limit t ^ oo. Our goal in this paper is to develop a 
the convolutionless generalized master equation suitable 
for a mesoscopic electron-phonon system using the fact 
that this equation becomes formally exact for arbitrary 
system sizes. In contrast to the derivation of the Pauli 
master equation, the only approximation that we will use 
is that of weak coupling. As a result we will arrive at a 
master equation with time dependent rate coefficients. 



where the prime at the summation sign indicates that the 
terms with n = m are excluded. This separation is use- 
ful for the following two reasons. First, in many systems 
only off-diagonal coefficients gnmi are small compared to 
gnni- Hence, V can be treated as a perturbation. Second, 
for many cases of interest, the initial density matrix com- 
mutes with Hq. In this case, the separation gives simpler 
forms of the master equations. 

For further analysis it is convenient to diagonalize Ho 
with respect to the normal mode degrees of freedom. 
This is achieved with the unitary operator 



U 
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as 



Hn = U- 
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where the renormaHzed electronic energies are 
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Applying the same unitary transformation to V gives 

X^' ( I \/ I t_L 
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two different master equations: the Nakajima-Zwanzig 
(NZ) equation ^JM|^ and the time-convolutionless (CL) 
master equation 



V = u-^vu 



(7) 



dVp 

dt 

dVp 
dt 



dTlC^''{t~T)Vp{T), (12) 



dr/C^^(r)Pp(i). (13) 



Introducing operators 



Mnmi = g-nmi « ■ + ^ 



,(8) 



The explicit expressions for superoperators JC^^ {t) and 
K'~^^{t) can be found in Ref. 3]. 

Even though Eqs. ^1 and are formally exact, it 
is impossible to calculate JC^^{t) and KP^{t) for most 
realistic systems. We will use the perturbation expan- 
sion of these superoperators. To the second order in the 
coupling constants 



we can rewrite 1^ in a compact form as 
= ^ \n){m\Mnmi- 



(9) 



The similar unitary transformations were applied to the 
electron-phonon Hamiltonians before [3i As can be 
seen from Eq. m the transformed picture the electronic 
transitions from state |n) to |m) are accompanied not 
only by the creation or annihilation of a single phonon of 
mode i but also by the displacements of all the normal 
modes. 

Let us take for the initial density matrix 



p{0) = \n){n\- 



Tr(e 



(10) 



where the electronic subsystem is prepared in a specific 
state \n) and the oscillators are in the displaced equi- 
librium configuration corresponding to this state. The 
transformed initial state p(0) = U^^pU has the form 



m = \n){n\ 



Tr(e~''S- 



(11) 



The electronic density matrix in the Eq. 1101 commutes 
with unitary operator Q , hence its time evolution will be 
the same in both the transformed and the untransformed 
pictures. For the purpose of developing a perturbative 
theory, it is more convenient at this point to work in 
with the transformed operators. 



III. MASTER EQUATION 

We now use the projection operator formalism to ob- 
tain the master equation for the diagonal part of the elec- 
tronic density matrix p*^' For a given projection 
operator V and the initial total density matrix that satis- 
fies p(0) = Vp{0), Vpit) can be shown to satisfy at least 



K.2 



(14) 



Here £o and Cv are the Liouville superoperators cor- 
responding to Hf) and V whose action on some density 
matrix p is given by 



Cop^ Hop- pHo, Cvp^Vp-pV. 



(15) 



In Ref. Q , Breuer and Petruccione show that to second 
order in the coupling constants, the convolutionless ex- 
peression (Eq. I13|l gives a better approximation to the 
exact solution than the Nakajima-Zwanzig equation (Eq. 
I12II . It also has an additional mathematical convenience 
of being local in time. Therefore, in the following analysis 
we will use the convolutionless approach. 

Since our initial state is given by Eq. ^Jwe will use the 
projection operator that acts on the total density matrix 
in the following way 



where 



Pp = 5]|n)(n|p:^Tr(|n)(n|p), 



(16) 
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Using the definition of /C2'^ (r) the following convolution- 
less equation is obtained for Pn = {'n\pei\n)^ 



dR 



Wnm{t)Pm - ^ W„„(t)P„ 



(18) 



The coefficients Wnmit) that are, in general, time depen- 
dent rates given by 



W™„(t) 
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and 



mnj 



-iff or 
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Due to the explicit form of operators M„,„j (Eq. IHJ the 
calculation of the correlation functions in Eq. [201 can be 
reduced to the averaging of the displacement operators 
over the equilibrium ensemble (Eq. I17|l . After straight- 
forward, but lengthy, calculations we obtain the principal 
result of this paper: 



X ((A„™(n, + l)e^'"'" 
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(27) 



In the case when all the diagonal electron/phonon terms 
vanish, gnni = 0, the correlation functions in Eq. 1221 re- 
duce to those obtained within the golden rule approach 
0. It is clear from Eq. [221 that for systems with 
gnni/'^i ^ 1 the goldcu rulc approach is not applicable. 

The qualitative behavior of the probabilities P„ {t) de- 
pends on the explicit form of the time dependent rates 
Wran{i) in Eq. 1181 Depending on a particular system, 
certain simplifications of Wmn{t) are possible. In par- 
ticular, for large systems we can expect the integrand in 



Eq. ^[to decay to zero and, correspondingly, Wmn{t) 
to become constant in the long time limit. If the typical 
electronic relaxation times are much longer then the de- 
cay rate of the integrand in Eq. ^1 then the Wmn{t) can 
be replaced by their time independent Markovian limits. 
Due to the rather complicated functional forms of qnm {t) 
and fnm{t) in Eq. 1221 we could not obtain an explicit ex- 
pression for Wrnn{t) in the Markovian limit. However, for 
models that satisfy the condition J^ji^j + h)^nmj ^ 
and with A„mi taken as smooth function of i, we can 
approximate fnm{T) in Eq. [^by a Gaussian, viz.. 



fnra ('^) ~ 6 



where 



3 



(28) 
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For uJiT 1, we can also approximate the sinw^r terms 
that appear in functions qnm{t) (Eq. I2fcif) by sinwiT w w^r 
so that 



At) 



where 



(30) 



(31) 



Under these approximations the long time (Markovian) 
limit of Wnra (i) Can be obtained as 



lim W„m(t) = y]i„ 



(32) 



where 



-i-ninjGnni{-(^i - ^3 + 7nm " + ^^ni) - {j^i + l)njGnm{(^i ~ ^3 + 7nm - + ^m) 
~{rij + l)niGmn{^j - + 7nm " + e^i)) 

-\-^nmi^7i'mjGnm('ynm ^m) 

)V (33) 
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FIG. 1: The population of the excited state as a function 
of time for the 80 mode system with Uj — Asin{j^), The 
sohd line corresponds to = 10, A = 0.1, p — —1/2, and 
e2 — ei — 1, the line of short dashes corresponds to ^ = 10, 
A — 0.2, p = 0, and £2 — ei = 1. The line of long dashes 
corresponds to the golden rule results for the Ohmic case. 
Temperature is 298 K for all cases. 

with 

Gmn{x) = J^^expl--^ — J. (34) 

IV. RELAXATION IN MODEL 
ELECTRON-PHONON SYSTEMS 

In what follows, we examine the solution of the Pauli 
equation for a number of model systems using the corre- 
lation functions (Eq. I22|) defined above. We take two ex- 
amples, one in consisting of a two state electronic system 
coupled to a single acoustic phonon branch with Ohmic 
and super-Ohmic couplings. The other taken from our 
recent work in describing secondary exciton formation at 
the domain boundary in organic polymer semiconductor 
heterojunctions. We also, then compare the results ob- 
tain using the present method to those obtained using 
Marcus' formulation of the rate-constant. 



A. Ohmic vs. Super-Ohmic Coupling 

As our first example we consider a system with two 
electronic levels coupled to N normal modes with fre- 
quencies uji = Asm{iTr/2N) with i G that model 
the acoustic phonon branch. We will choose the coupling 
coefficients in the following form gi2i — g2ii — giu — 
X(uji)P and (7221 = — A(wi)^, where A and p are parame- 
ters. We will consider the cases oi p = —1/2 and p — 
that correspond approximately to the Ohmic and super- 
Ohmic baths in the golden rule approach to relaxation. 

Taking the initial state to be prepared in the higher- 
lying electronic state, electronic population of the higher 
electronic state as a function of time is shown on Fig. ^ 



FIG. 2: Transition rates Wi2{t) and W2i{t) for the super- 
Ohmic model 



It can be seen that the relaxation does not follow a simple 
exponential decay. This behavior can be understood by 
looking at the explicit form of the rate coefficients Wi2{t) 
and W2i{t) as functions of time as shown on Fig. |21for 
the super-Ohmic case. Only after the rate coefficients be- 
come nearly constant in time does exponential relaxation 
occur. For the chosen parameter values, substantial non- 
exponential relaxation occurs before that time. Note the 
for the chosen parameters the complete electronic relax- 
ation occurs before the onset of any possible recurrence 
phenomena due to the finite number of modes. Such re- 
currences will occur at the time scale of the order of 2t: 
divided by the smallest normal mode frequency. It can 
also be seen on Fig. 2 that coefficients Wi2{t) and W2i{t) 
can acquire negative values. This is in contrast to the 
golden rule approach where the coefficients are always 
positive and time independent. Because of the time de- 
pendence of coefficients Wi2{t) and W2i{t) and their ap- 
proximate nature the time-convolutionless master equa- 
tion is not guaranteed to preserve the positivity of prob- 
abilities. Fig. 2 also shows the electronic relaxation for 
the Ohmic model obtained by applying the golden rule to 
the Hamiltonian |^ when coefficients gnni are treated on 
the same basis as gnmi {n 7^ ^n). No reliable golden rule 
results can be obtained for the super-Ohmic case because 
of the relatively small number of modes in the model. 



B. Exciton regeneration at a semiconductor 
interface 

As another example we consider the electronic relax- 
ation in the conjugated polymer heterojunctions previ- 
ously investigated in Ref. [lOl llll by a different ap- 
proach. Here we consider the decay of an excitonic state 
into a charge-separated state at the interface between 
two semiconducting polymer phases (TFB) and (F8BT). 
Such materials have been extensively studied for their po- 
tential in organic light-emitting diodes and organic pho- 
tovoltaics [la. IT3 . [l5l | . At the phase boundary, the ma- 
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FIG. 3: The population of the excited state as a function 
of time for the TFB:F8BT heterojunction for three differ- 
ent temperatures obtained by solving the time-convolutionless 
master equation (ME) or through the Marcus-type approxi- 
mation (MT). 

terial forms a type-II semiconductor heterojunction with 
the ofF-set between the valence bands of the two materials 
being only slightly more than the binding energy of an 
exciton placed on either the TFB or F8BT polymer. As 
a result, an exciton on the F8BT side will dissociate to 
form a charge-separated (exciplex) state at the interface. 

F8BT* : TFB — > F8BT" : TFB+. 

Ordinarily, such type II systems are best suited for photo- 
voltaic rather than LED applications However, LEDs fab- 
ricated from phase-segregated 50:50 blends of TFB:F8BT 
give remarkably efficient electroluminescence efficiency 
due to secondary exciton formation due the back-reaction 

F8BT" ; TFB+ — > F8BT* : TFB. 

as thermal equilbrium between the excitonic and charge- 
transfer states is established. This is evidenced by long- 
time emission, blue-shifted relative to the emission from 
the exciplex, accounting for nearly 90% of the integrated 
photo-emission. 

Here, as above, as consider only two lowest electronic 
levels corresponding to 

\XT) = F8BT* : TFB & \CT) = FSBT' : TFB+. 

As was shown in Ref. 0, IT^ IT^ for such systems 
there are two groups of phonon modes that are coupled 
strongly to the electronic degrees of freedom as evidenced 
by their presence as vibronic features in the vibronic 
emission spectra, namely: low frequency torsional modes 
along the backbones of the polymer and higher frequency 
C=C stretching modes. The normal modes within each 
group have about the same frequencies and form two sep- 
arated bands. The interested reader is referred to Ref. 
[T^ IT^ IT^ for specific details regarding the parameteri- 
zation of the model. 



FIG. 4: Transition rates Wi2{t) and W2i{t) for the 
TFB:F8BT heterojunction (see Ref. for the model de- 

tails) involving 24 normal modes at 298 K. 



As before, we consider the initial state as being pre- 
pared in the \XT) state corresponding to photoexcitation 
of the F8BT polymer. The results for the higher elec- 
tronic level population obtained by numerically solving 
Eq. H18|l are shown on Fig. 13 The time dependence of 
coefficients Wi2{t) and W2i(t) is shown on Fig. ^ As in 
the previous example the electronic relaxation does not 
follow a simple exponential pattern. The main difference 
between this model and the previous example is that the 
electronic relaxation time is of the same order as the re- 
currence time for coefficients Wi2{t) and W2i{t). As can 
be seen from Fig. ^ the time dependence of these coef- 
ficients has the form of approximately constant regions 
abruptly changed at regular recurrence intervals. This 
type of time dependence makes the electronic relaxation 
look like a series of exponential relaxations with changing 
rates. As in the previous example, coefficients Wi2{t) and 
W2i{t) can become negative as the relaxation proceeds 
towards the equilibrium population. 

Note also, that the because the relaxation does not 
obey a simple exponential rate law, the initial decay 
slightly "overshoots" the final equilibrium population. 
This is most evident in the highest temperature case 
considered here (T = 340K). Since the \XT) is also the 
emissive state, photo-emission (not included herein) de- 
pletes this population on a nanosecond time scale (radia- 
tive lifetime). Population back-transfer from the \CT) 
to maintain a thermal equilibrium population then leads 
to the continuous replenishment of the emissive species 
such that nearly all of the CT population contributes to 
the formation of secondary (regenerated) excitonic states 

naiiiiiiiiiiiiiii. 

It is of interest to compare the relaxation dynamics in 
TFB:F8BT heterojunction obtained through the appli- 
cation of the time-convolutionless master equation with 
other approaches. As an example we will consider the 
Marcus-type approach which is widely used to study elec- 
tron transfer in chemical systems 0, 0, |20| . Be- 
cause of the coordinate dependency of the coupling, some 
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care must be taken in deriving the Marcus rates and 
we give details of this in the appendix. The results 
for the relaxation dynamics using these rates applied to 
TFB:F8BT heterojunction for three different tempera- 
tures are shown in Fig. |3| It can be seen that apart from 
a more complicated time dependence the master equation 
approach gives faster relaxation when compared to the 
Marcus-type picture. The discrepancy between the two 
approaches can be explained by the fact that the Marcus 
approximation is assumed to be valid when kT ^ hoji for 
all the normal modes. In the case of the TFB:F8BT het- 
erojunction this condition is not satisfied for the higher 
frequency modes. 



foundation. Seed funds from the Texas Center for Su- 
perconductivity (TcSUH) are also acknowledged. 



APPENDIX: DERIVATION OF MARCUS- TYPE 
RATE EXPRESSIONS 

To apply the Marcus approach for the system at hand, 
it is convenient to rewrite Hamiltonians l(2Jl and |(2Jl in 
terms of the coordinate qi and momentum pi operators 
for the normal modes as 



V. CONCLUSIONS 

We derived the explicit form of the time- 
convolutionless master equation for the electronic 
populations in the electron-phonon systems. We applied 
it to study the electronic relaxation in the two-level elec- 
tronic systems coupled to the Ohmic and super-Ohmic 
baths, as well as in the realistic conjugated polymer 
heterojunction. The time evolution of the electronic 
populations is, in general, more complex then the one 
obtained with the Pauli master equation and depends 
on the explicit form of the time dependent coefficient. 
The time-convolutionless master equation can account 
for the appearance of recurrence effects in smaller 
electron-phonon systems. 
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n ni i i 

and an off-diagonal part V 

V = ^ 5„m^|"-)("^|g^, (A.2) 
nmi 

where 

9nmi — 'V'^^idnmi- {-^■^) 

The coordinate and momentum operators are then 
treated as static parameters and the rate constants for 
the transition from state \n) to state |m) is given by av- 
eraging the golden rule transition rates over the initial 
equilibrium distribution of the coordinates corresponding 
to the equilibrium ensemble for Hamiltonian Hq assum- 
ing the that the electronic state of the system is \n) 
i.e., 



d{q^}f{{q^})\Vnmi{q^})\^S{Un{{q^}) " U„,{{q,}) 



(A.4) 



where {qi} denotes all coordinate variables and 



(A.5) 



and 



Vnmiiqi}) = '^gnmiqi- (A.6) 

i 

Here Z is the partition function and Un is given by 

£n ^~ / , 9nniQi 



Note that Eq. I A. 41 differs from the similar expression of 
Ref. by the presence of coordinate dependence in the 
interaction parameters Vnm{{qi})- The integrations in 
Eq. I A. 41 can be explicitly performed giving 



/ ^^nraPi 



nnri^ run 



\ nT7 



I TT/3 



■ exp 



(A., 
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Here AEnm and E^„^ are given by 



■i * 2 * 



(A.IO) 



and are usually referred to as the driving force and reor- 
ganization energy, respectively. Parameters Pnm, Fnm, 
and Hnrri are given by 

p QnmiiQnni 9mmi) /a \ 



Hn 



2u;f 



jnrni 



(A.12) 
(A.13) 



It can be verified from Eqs. (|A.4p or l|A.8|l that rates k„ 
satisfy the principle of the detailed balance, i. e.. 



k. 



h 



(A. 14) 
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